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We present the first Bayesian constraints on the single field inflationary reheating era obtained 
from Cosmic Microwave Background (CMB) data. After demonstrating that this epoch can be fully 
characterized by the so-called reheating parameter, we show that it is constrained by the seven 
years Wilkinson Microwave Anisotropies Probe (WMAP7) data for all large and small field models. 
An interesting feature of our approach is that it yields lower bounds on the reheating temperature 
which can be combined with the upper bounds associated with gravitinos production. For large 
field models, we find the energy scale of reheating to be higher than those probed at the Large 
Hadron Collider, > 17.3 TeV at 95% of confidence. For small field models, we obtain the two- 

1/4 

sigma lower limits p T £ > 890 TeV for a mean equation of state during reheating u> rc h = —0.3 and 
Pteh > ^90 GeV for = —0.2. The physical origin of these constraints is pedagogically explained 
by means of the slow-roll approximation. Finally, when marginalizing over all possible reheating 
history, the WMAP7 data push massive inflation under pressure (p < 2.2 at 95% of confidence 
where p is the power index of the large field potentials) while they slightly favor super-Planckian 
field expectation values in the small field models. 

PACS numbers: 98.80.Cq, 98.70. Vc 



I. INTRODUCTION 



The current ongoing flow of high accuracy astrophysi- 
cal observations has important consequences for our un- 
derstanding of the very early Universe. Inparticular, the 
widely accepted inflationary paradigm [l|-[4[ (for a review, 



see e.g. Refs. [5-8]) is now under close scrutiny. Accord- 
ing to this scenario, the Cosmic Microwave Background 
(CMB) anisotropies and the large scale structures orig- 
inate from the unavoidable quantum fluctuations of the 
inflaton and gravitational fields in the very early Universe 
subsequently amplified during inflation |9l-fl3j. One can 
show that the corresponding power spectrum of the cos- 
mological fluctuations naturally acquires an almost scale 
invariant form which is fully consistent with all observa- 
tions. Another crucial property of the inflationary power 
spectrum is that the slight deviations from scale invari- 
ance are linked to the microphysics of inflation [T3 - IT^ . 
Therefore, by measuring these deviations, one can probe 
the shape of the inflaton potential and, therefore, learn 
about the physical origin of the inflaton field. 

It is often claimed from the above properties that ob- 
servations give access to a limited part of the potential 
only, namely the one which is slow-rolled over by the 
inflaton when scales of astrophysical interest today left 
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the Hubble radius. This observational window represents 
a range of approximately 7 e-folds or three decades in 
wavenumbers. However, inflation does not consist of the 
slow- roll phase only and the pre and/or reheating period 
is also of fundamental importance since it allows us to un- 
derstand how inflation is connected to the hot big-bang 
phase (l7l - i20j . This physical phenomenon is related to 
a different part of the inflationary potential, usually the 
one located close to its true minimum, i.e. a few decades 
in e-folds away from the observable window. 

Observation of pre- /reheating effects can be achieved 
in two ways. First, the power spectrum can evolve on 
large scales when the inflaton field oscillates around the 
minimum of its potential. However, this happens only 
in quite complicated models, typically those containing 
more than one field [2l| - i23j . In fact, it was recently 
shown that this type of effect can also happen in single 
field inflation but on much smaller scales [24-26]. Sec- 
ond, the duration of the pre- /reheating phase can signif- 
icantly modify the position of the observational window 
mentioned above. Put differently, at fixed astrophysi- 
cal scales today, changing the pre-/reheating duration 
is equivalent to moving the window along the poten- 
tial, hence probing different values of the power spec- 
trum spectral index, amplitude of the anisotropies and 
tensor-to-scalar ratio. Obviously, this cannot be done 
arbitrarily because CMB data impose accurate bounds 
on their value. Conversely, this opens up the possibility 
to constrain the pre- /reheating duration and/or its equa- 
tion of state from CMB data [23. Notice that a direct 
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detection of primordial gravitational waves would also al- 
low us to probe the reheating temperature, as shown in 
Refs. [H,[29|. 

The goal of this article is to address this question for 
the standard scenarios of inflation. It is traditional to 
study three categories of models usually considered as 
representative of the full inflationary space. These mod- 
els are large field [3(| , small field @, Q , and hybrid in- 
flation [3l| . Hybrid scenarios involve multiple fields and, 
therefore, the power spectrum can change during the pre- 
heating phase. This makes this class of scenarios deserv- 
ing of a separate investigation. For this reason, in this 
article, we limit ourselves to the class of large and small 
field models. 

In the following, we will use the term "reheating" to 
refer to the pre-/reheating phases of the Universe de- 
fined to have occurred just after the end of inflation and 
just before the radiation dominated era. So far, the con- 
straints on the reheating energy scale are not so numer- 
ous. Obviously, it should be less than the energy scale 
of inflation which implies that T re h < 10 16 GeV. In ad- 
dition, if one assumes that supersymmetry is the cor- 
rect extension of the standard model of particle physics, 
then constraints from Big-Bang Nucleosynthesis (BBN) 
on unstable gravitinos lead to a reheating temperature 
Troii < 10 7 GeV HHHj]. Notice that this constraint 
can nevertheless be avoided if one considers the scenario 
of Ref. [il] . Reheating itself should also proceed before 
BBN and this implies that T reh > 10 MeV. We see that 
the reheating temperature is poorly constrained, in par- 
ticular, its lower limit. As a matter of fact, the work pre- 
sented here precisely yields a lower limit on the reheating 
energy scale from the current seven years Wilkinson Mi- 
crowave Anisotropies Probe (WMAP7) data [4^ - [46j . 

In order to derive constraints on the reheating phase, 
we make use of Bayesian techniques and utilize a full 
numerical approach j42j. This has several advantages. 
First, it is exact and rests only on the linear theory of 
cosmological perturbations: the method remains accu- 
rate when the slow-roll approximation breaks down, as 
one expects near the end of inflation. Second, and of 
particular importance for the present work, it permits a 
new treatment of reheating. Indeed, instead of viewing 
the reheating parameters as nuisance parameters, they 
can easily be included in the Bayesian data analysis pro- 
cess. Third, the evolution of cosmological perturbations 
in the hot big-bang eras already relies on numerical codes. 
Treating perturbations during inflation in the same way 
allows the whole procedure to be automatized and to be 
easily extended to other scenarios. Fourth, the numer- 
ical approach allows us to address the question of the 
priors choice in a particularly well-defined way. Indeed, 
from a physical point of view, our prior knowledge is on 
the inflationary theory and not on the shape of the pri- 
mordial power spectra which is actually a model predic- 
tion. Therefore, it is better, and easier, to choose prior 
probability distributions directly on the model parame- 
ters, such as the power index of the large field potentials. 



This reflects the fact that a model of inflation is not a 
disembodied mathematical structure that one only needs 
to "fit" but a physical scenario rooted in high energy 
physics that one needs to understand. 

This paper is organized as follows. In Sec. [TT1 we ex- 
tend the above discussion and explain in detail why the 
reheating epoch can be constrained with CMB data. In 
particular, we introduce the so-called reheating param- 
eter which depends on the reheating duration and on 
the mean equation of state of the fluid dominating the 
Universe during this epoch. Then, using the slow-roll 
approximation, we analytically demonstrate that the ac- 
curacy of the WMAP7 data is now sufficient to obtain 
some constraints on the reheating era. In Sec. IIII1 us- 
ing a full numerical integration of the tensor and scalar 
power spectra coupled to Bayesian methods, we derive 
the constraints that any reheating model has to satisfy. 
Then, assuming specific values for the mean equation of 
state, we translate these constraints into new lower limits 
for the reheating energy density and/or reheating tem- 
perature. These results significantly improve the bounds 
coming from the Big-Bang Nucleosynthesis. In Sec. IIV1 
we recap our main findings and discuss how our results 
are modified by the inclusion of others CMB data sets. 
In Appendix [A[ we work out a typical example which il- 
lustrates the robustness of our assumptions: a noninstan- 
taneous transition between reheating and the radiation 
dominated era when one considers the finite decay width 
of the inflaton field. Finally, as a by-product of our data 
analysis, Appendix iBlpresents the updated WMAP7 con- 
straints on the spectral index, tensor-to-scalar ratio and 
first order slow-roll parameters marginalized over second 
order effects. 



II. PHYSICAL ORIGIN OF THE CONSTRAINT 

Before presenting and discussing the constraints on the 
reheating temperature, we explain why and how these 
ones can be inferred from high accuracy CMB observa- 
tions. In particular, we use the slow-roll approximation 
to explicitly illustrate the method. 



A. Parametrizing the reheating 

The evolution of scalar (density) perturbations is con- 
trolled by the so-called Mukhanov-Sasaki variable v k . If 
matter is described by a scalar field (as is the case during 
inflation and pre-/reheating), then its equation of motion 
is given, in Fourier space, by @-H, HU 



v k = 0. 



(1) 



Here, a prime denotes a derivative with respect to con- 
formal time. The quantity k is the comoving wave num- 
ber and ei = —H/H 2 is the first Hubble flow func- 
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tion [4g], H = a/ a being the Hubble parameter and 
a the Friedmann-Lemaitre-Robertson-Walker (FLRW) 
scale factor (a dot means derivative with respect to cos- 
mic time). The quantity Vk is related to the curvature 
perturbation Oc through the following expression: 



1 



(2) 



where M pl stands for the reduced Planck mass. As a 
consequence, the power spectrum of can be expressed 
as 



P c (k) 



fc3 \c 1 2 
— 2 ICfcl 



fc 3 



2tt 



4ir 2 Ml 



1% 



(3) 



In order to calculate V((k), one needs to integrate 
Eq. (TTJ), which requires the knowledge of the initial con- 
ditions for the mode function Vk- Since, at the begin- 
ning of inflation, all the modes of astrophysical interest 
today were much smaller than the Hubble radius, the ini- 
tial conditions are chosen to be the Bunch-Davis vacuum 
which amounts to 



lim Vk 



1 



-ihr\ 



2k 



(4) 



where r\ denotes conformal time and % = aH is the con- 
formal Hubble parameter. The importance of the curva- 
ture perturbation lies in the fact that it is directly related 
to CMB anisotropics, the two point correlation function 
of which can be expressed in term of the spectrum of Qk ■ 
Moreover, under very general conditions (including the 
assumption that inflation proceeds with only one field), 
Cfc is a conserved quantity on large scales and, there- 
fore, can be used to propagate the inflationary spectrum 
from the end of inflation to the post-inflationary era [4i| . 
In other words, the power spectrum is not affected by 
the post-inflationary evolution, in particular by the pre- 
/reheating epoch. 

However, this does not mean that the reheating era 
has no effect on the inflationary predictions. On the con- 
trary, the relation between the physical scales at present 
time and during inflation depends on the properties of 
this phase of evolution. As a consequence, in order to 
calibrate the inflationary spectrum with respect to the 
physical scales of astrophysical interest today, it is neces- 
sary to know how the reheating phase proceeded. Con- 
versely, this also opens the possibility to constrain the 
physical conditions that prevailed at that time by means 
of CMB observations. 

In order to put the above considerations on a quantita- 
tive footing, let us rewrite Eq. ([1]) in terms of the number 
of e-folds during inflation, N = ln(a/ai n i), where ai nl is 
the value of the scale factor at the beginning of inflation. 
It takes the form 



d 2 v k 
dN 2 



1 dH dv k 
% dN dN 



U S (N) 



v k = 0, (5) 



where U S (N) is an effective potential for the perturba- 
tions which depends on the scale factor and its deriva- 
tives only. All the terms in this equation but k/T-L are 
completely specified by the inflationary background evo- 
lution. In practice, we are given a physical scale today, 
say k/a now (for instance k/a now — 0.05 Mpc -1 ) and we 
need to express k/T-L in terms of k/a now and quantities 
defined during inflation. Straightforward considerations 
lead to 



k 

n 



H(N) 



(0) 



where N T is the total number of e-folds during inflation 
and Tfc is defined by 



(1 + Zend) , 



(7) 



with z cn( j being the redshift of the end of inflation. As 
expected depends on the whole post-inflationary his- 
tory through Zend- During this post-inflationary history, 
only the reheating phase is poorly known and represents, 
by far, the main source of uncertainty for the inflationary 
predictions. For convenience, we rewrite as 



T fc = 



Pcnd 



1/4 



R 



rad ' 



(8) 



thus defining the new parameter i? ra( j. This parame- 
ters plays a crucial role in this article. In the above 
equation, p cn( j is the energy density at the end of in- 
flation, p cr ; is the present day critical energy density and 
f2 7 ~ 2.471 x 10 -5 /i -2 is the density parameter of ra- 
diation today. As a result f2 7 p C ri = P7 is the present 
day radiation energy density and does not depend on h 2 . 
The above equations make clear that the parameter i? ra d 
must be specified if one wants to compare an inflationary 
model to observations. 

In fact, the quantity i? ra( j has a simple physical in- 
terpretation. Let us assume that the reheating phase 
is dominated by a conserved effective fluid with energy 
density p and pressure P. The fact that we assume the 
effective fluid to be conserved is not a limitation. For in- 
stance, in a simple model where the inflaton scalar field is 
coupled to radiation (see Appendix [A"|l . the effective fluid 
is just defined by p — + p 7 and P = + p 7 /3. The 
scalar field and the radiation are not separately conserved 
but the effective fluid is. Then, it is straightforward to 
show that 



p(N) = p ond exp |-3 J [1 + w rch (n)] dnj 



(9) 



where U! rcn = P/p is the equation of state function during 
reheating. Using this expression, one obtains 



AiV 

hli? rad = (-1 + 3lO re h) 



(10) 
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where 



where we have defined 



AN = N rch - N T 



(11) 



is the total number of e- folds during reheating, -ZV ro h being 
the number of e- folds at which reheating is completed and 
the radiation dominated era begins. The quantity w Te h 
stands for the mean equation of state parameter 



1 



W Ie h 



AN 



Nr. 



— — / w leh {n)dn. 



(12) 



A? 



Therefore, the parameter i? ra a only depends on what 
happens during reheating. To put differently, it singles 
out in the expression of the contribution coming from 
reheating. Equation (|10[) also allows us to understand 
why i? rad carries the subscript "rad". Indeed, if the ef- 
fective fluid is equivalent to radiation, then w rc h = 1/3 
and In i? ra d = 0. The physical interpretation is very clear: 
in this case the reheating stage cannot be distinguished 
from the subsequent radiation dominated era and, there- 
fore, cannot affect the inflationary predictions: as a con- 
sequence i? ra d = 1 in Eq. ((SJ. 

In fact, one can even go further and express i? ra d in an 
even more compact form. Using Eq. ([5]), one can write 
Proh = Pcnd exp [— 3AiV(l + Wrch)] from which, together 
with Eq. (fT0|). one obtains 



InR 



rad 



1 - 3w rG h 
12(1 + W ieh ) 



In 



Prch 
Pond 



(13) 



where p ren has to be understood as the energy density at 
the end of the reheating era, i.e. p(iV r ch)- 

Let us summarize our discussion. In order to calculate 
the power spectrum of the inflationary cosmological per- 
turbations, one needs to solve Eq. ([5]). In this formula, all 
the terms are accurately known during inflation except 



k 
H 



Pend 
Pi 



1/4 



1 



-,N„ —N 



H(N)R r , 



ad 



(14) 



and the theoretical uncertainty in this expression solely 
comes from the parameter i? ra d which depends on reheat- 
ing only (more precisely, on the energy density at the end 
of reheating, /3 rcn , and the mean equation of state w ro h)- 



B. Why CMB observations constrain reheating 

Having discussed the physical interpretation of i? ra d, 
we now explain how the CMB observations can constrain 
its value. For this purpose, we reexpress i? ra d in terms 
of quantities defined at the Hubble radius crossing. One 
obtains 



lni? rad = N T — AT* + Nq — - In 



1, ( H * 



-m 

4 



3 V^cnd 3 — ei* 

£l* U* 3 — £l end 



(15) 



N = In 



fc/dn 



1/4 

p 7 



(16) 



In this formula, N* is the e-folds number at which the 
scale k/a now crossed out the Hubble radius during infla- 
tion (all the quantities with a subscript "*" are evaluated 
at that time) and V(cf>) is the inflaton potential. Despite 
the appearance of the first Hubble flow function, this 
equation is exact (moreover, we also have ei on d = 1). At 
leading order, one has 



Ml 



c 1* 



8tt 2 P*, 



(17) 



where the amplitude of the scalar power spectrum at the 
pivot scale P* = V$ (A;*) is directly related to the Cosmic 
Background Observer (COBE) normalization. 

The above equation (1151) can be used in two different 
manners. The first way is to assume something about 
i? ra d and to derive the corresponding range of variations 
of the inflationary slow-roll predictions N* and ej(iV*). 
In other words, this determines how the inflationary pre- 
dictions depend on the details of the reheating era. This 
approach is the one usually considered in the literature 
to compare inflationary predictions to the current con- 
straints on the slow-roll parameters e;* (or spectral index 
and tensor-to-scalar ratio). Unfortunately, the assump- 
tions on i? ra d are rarely explicit and comparison is only 
made by choosing reasonably assumed values of iV» : typ- 
ically 30 and 60 e-folds as one may derive under generic 
assumptions However, as Eq. (fT5")) explicitly shows, 
once V(4>) is chosen, and the tilt and amplitude of the 
scalar perturbations measured, TV* is directly related to 
Pradj which itself, as already noticed, depends on the en- 
ergy density p ren at which reheating ends and z/> rcn . As a 
result, the range of variation for N* can only be known 
once a reheating model is assumed. Without such an as- 
sumption, from one model to another, an assumed value 
of N„ may inconsistently imply that the reheating occurs 
after nucleosynthesis, or even at energy densities higher 
than Pend' This type of model would therefore appears to 
be compatible with the CMB data favored power spectra 
while being totally inconsistent with standard cosmology. 

Let us now see how it works in practice. In order to be 
consistent with the standard cosmological model, In i? ra d 
cannot take arbitrary values. One should have u> rcn < 
1 to respect the positivity energy conditions of General 
Relativity and w rc h > — 1/3 by the very definition of 
reheating which is not inflation. Notice that we impose 
conditions on the mean value of the equation of state 
only. In addition, reheating should occur after inflation 
and before BBN, i.e. p nuc < p re h < p cn d, with 



(lOMeV) 4 



(18) 



This allows us to explicitly use Eq. (fT5|) . Combined 
with Eq. (|13p . we can determine the range of variation 
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of AJV, = N T - N, G [AiVJ 
manipulations lead to 



A7V™ d ]. Straightforward 



AN" 



-N + In 



In 



Pcnd 



M Pi ; 3(l + iD rch ) M* 



1 - 3w rch . /Cnuc 

In ■ 



12(l + tD rcn ) M* ' 
while if one chooses p rc h = Pcnd, one obtains 

Aiv cnd = _ iVo + ln f _ 1 ln Pcnd 



(19) 



\M r J 4 M* 



(20) 



Interestingly enough, the last equation no longer depends 
on w) re h- This is of course because requiring p ro h = p e nd 
means that one immediately reheats the Universe after 
inflation. It is important to notice that these equations 
are algebraic for A7V" UC and A./V° nd because £T* and p e nd 
are also functions of AiV*. The corresponding range of 
variations of the inflationary predictions is determined 
by calculating Ae« = ej(AiV*) with AA* given above. 
To proceed further, one needs to specify the model of 
inflation. In the next section, one considers the proto- 
typical scenario of chaotic inflation as well as the small 
field models. 

Before dealing with these explicit examples, let us 
briefly anticipate and discuss the second way of using 
Eq. (TT5|) . It consists of considering i? ra d as an observable 
model parameter and in including it in the data analysis, 
as it should be from a Bayesian point of view. If we are 
given a specific potential, then V on d is explicitly known. 
CMB data put a limit on H%/eu through the amplitude 
of the anisotropies, as well as on ei* from the tensor-to- 
scalar ratio. As a result, one expects CMB data to also 
give some information on i? ra d- This is the subject of 
Sec. IIIII in which we perform a Bayesian data analysis 
of the WMAP7 data for both the large and small field 
models by including the reheating. For the first time, 
we find that -R ra d is not only a nuisance parameter for 
inflation but ends up being constrained by the WMAP7 
data. We then discuss the physical implications of these 
bounds and show that CMB data give us a lower bound 
on the energy scale at which the reheating ended. 



C. Large field models 

We now consider the archetypal model of inflation, 
namely, large field inflation. This working example is 
important because it allows us to show explicitly which 
type of constraints one should expect. Large field models 
are characterized by the potential 



V((f>) = M 4 



(21) 



where M is an energy scale which fixes the amplitude of 
the CMB anisotropies and p is a free index. In this case 



the slow-roll trajectory is explicitly known and one can 
calculate </!>* , the field vacuum expectation value at Hub- 
ble radius crossing from (fiend /M p[ — p/y/2, the field vac- 
uum expectation value (VEV) at which inflation stops. 
One gets [27lj 



= 2pM*AJV. 



^end- 



(22) 



The reheating phase in large field models proceeds by 
parametric oscillations around the minimum of the po- 
tential and it is well known that the corresponding equa- 
tion of state parameter is given by @, [13, E, IHfjj 



P 



P- 



(23) 



In particular, for p = 2, one obtains w) re h = 0: that is to 
say the oscillatory phase is equivalent to a matter domi- 
nated era (the quartic case corresponding to a radiation 
dominated era, and so on). Although this formula is de- 
rived without taking into account the coupling between 
the inflaton field and radiation, we show in Appendix lAl 
that it is a very good approximation. 

Knowing explicitly the equation during reheating, we 
are now in a position where the algebraic equations (|19p 
and (j2"0|) can be solved exactly. After some algebra, one 
obtains 



AA" UC = - l - 



p p 2 -2p + 4 



x exp 



Up 
12p (AA n 



W 



Up 



P 
P/4)' 



2p + 4 



p 2 - 2p + 4 



(24) 



and 



p p 



x exp 



-Wo4 — 



8 u [ p-2 

8 (A/" cnd + p/A) 
P~2 



(25) 



where Wo is a Lambert function. Both quantities J\f nuc 
and J\f cnd depend only on the model parameter p and the 
amplitude of the observed anisotropies. Explicitly, J\f nuc 
reads 




where the amplitude of the CMB anisotropies has been 
expressed in terms of the quadrupole moment 



Q 



rms-PS 



5C 2 

47T 



6 x 10 



-6 



(27) 



In Eq. (|26p . the last term vanishes for p — 4 since, as 
already noticed above, the phase of oscillations is equiv- 
alent to a radiation dominated era which cannot be dis- 
tinguished from the subsequent hot big-bang epoch. On 
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the other hand, the constant AT ond can be expressed as 



# cnd = -N + - In ^27rVl20 
1 , / 2 x -v 



:Qrms- 



PS 



T 



(28) 



From Eqs. (p4|) and (|25|) , we immediately deduce that, 
for the large field models, the range of allowed val- 
ues of AiV* strongly depends on p. In Figs. Q] and [5J 
we have plotted the large field predictions, obtained 
from Eqs. ([24]) and (|25|) . for the slow-roll parameters 
ei(AN^) and £2 (A A*) compared to the one- and two- 
sigma WMAP7 confidence intervals (see Appendix IB"]) . 
The first figure represents the slow-roll predictions in the 
plane (n s ,r), while the second one corresponds to the 
plane (61,62). The annotated values trace the quantity 

l°g(9*^ 4 ^rch/GeV), where the reheating temperature is 
defined by the relation 



i/V f 30 



1/4 



(29) 



and g* is the number of relativistic degrees of freedom 
at that time. Large values of p cannot explain the cur- 
rent measurements of n s and r, while the low values can. 
Therefore, it is clear that there now exists a lower bound 
on the reheating temperature. In Sec. IIII1 we derive the 
Bayesian two-sigma limits on p ro h (or gV^T^) by in- 
cluding the reheating parameter into the data analysis 
process. These plots should make evident that the re- 
heating in the large field models is already observable 
with the current CMB data, and more than being a nui- 
sance parameter, it is actually constrained. 

Let us also remark that the case p — 4 is particularly 
interesting; see the blue point annotated "16" in Figs. Q] 
and [2] Indeed, the value p = 4 is the extreme case in 
which is actually fixed to 



ANP= 4 = -N + lln 



rms— PS 



58.5. (30) 



This is why this model is represented by a single point 
in Figs. [T] and [2] Making any other choice is equivalent 
to assuming a more complicated reheating model which 
should at least be specified. For instance, Ref. (see 
Fig. 19) uses two values, 50 and 60 e-folds, instead of one. 
From the above considerations, it is clear that 50 is much 
too small. But, of course, one can always assume that 
the shape of the potential in the slow-roll regime is not 
the same as in the reheating regime (actually, this has 
to be the case for small field models, see below). In this 
case V((f>) oc 4 is not relevant during the oscillations of 
the field and w rc h 7^ 1/3. However, as discussed in the 
next section, even in this case, the reheating epoch is 
still constrained. The same range of variations for AA* 
has also been used in Ref. [5l[ (see Fig. 2) for the <fi 2 
model. Compared to our results, 50 is too high and ex- 
cludes models which are still allowed while 60 predicts a 



reheating energy scale higher the energy scale at the end 
of inflation: p reh > p end . 



D. Small field models 

One of the reasons leading to such a strong reheating 
influence on the large field model predictions comes from 
Eq. (|23[) . Once the potential is chosen, the spectral index 
and tensor-to-scalar ratio are intimately linked to the way 
the reheating proceeds. One may therefore wonder how 
the reheating can influence the model predictions in a 
case where it is unrelated to the shape of the primordial 
power spectra. As a motivated example, we discuss in 
this section the case of the small field models ending with 
a reheating characterized by a mean equation of state 
u) ren . The small field potential reads 



V{4>) = M 4 



(31) 



where /1 represents a VEV for the field cf>, p a power index, 
and M fixes the amplitude of the observed anisotropics. 
Inflation proceeds from small to large values of the field. 
For convenience, we denote by x the field value in units 
of jit, i.e. x = </>//•*• The slow-roll trajectory in terms of x 
reads 123 



AN^^lxl 



2pM* t 



2 2-p _ 2 1 ^ 2 
X* Xcnd ' 



p-2 



2-p 



Xcnd 



where, again, Xcnd is defined by ei(xcnd) = 1, i-e. 



V2 



p M pl 



Xcnd) ' 



(32) 



(33) 



Both of these equations do not have an explicit solution 
(unless p = 2) but can be numerically solved for a given 
set of model parameters fj, and p. For this reason, instead 
of deriving the reheating allowed values for AiV*, it is 
more convenient to derive the reheating allowed values 
for x*- 111 iac t> X* should lie between x* uc an d X* nd , 
the field values such that reheating ends, respectively, at 
BBN and just after inflation. After some algebra, one 
finds x* uc to be the solution of 



X* 



2-p 



2pMl \- 'p-2 

_ 3w re h + 1 

3 + 3w re h 
,2 



2pM 2 



M 'x 2 , 

Acnd 



p-2 



3w re h 

3 + 3w rch 



2-p 
Xcnd 



J- nuc , 



(34) 
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FIG. 1: Reheating consistent slow-roll predictions for the large field models in the plane (n s , r). The two contours are the one 
and two-sigma WMAP confidence intervals (marginalized over second order slow- roll). The two lines represent the locus of the 
p > 1 and p = 2 models while the blue point annotated "16" corresponds to p = 4. The annotations trace the energy scale at 
which the large field reheating ends and correspond to log(<?y 4 T re h/GeV). Clearly, these values are limited from below to stay 
inside the two-sigma contours. 



with 



jt-nuc = _ N 



1 + 3w rch 




3 + 3w rc h 



1 _ o, T /n 1 / 4 ' 
i/i p \ , 1 •J w rch , / Pnuc 

111 (1 - Xcnd) + — ln 



3 + 3w rc h 



(35) 

Similarly, solving for p rc h = p en( j gives xt nd as the solu- 
tion of 



Xl + ~^xl~ P ) + \ ln(l ~xl)-\ Hxl~ l ) 
p — 2 ) 4 2 



2pM 2 



Xend 



p-2 



X c 2 n 7 +^ nd 



(36) 



where 



-JVb + \ In f2^7T20 Qrn "- PS 

Z \ 1 



iln 
4 



(37) 



As a result, for given values of fi, p and w ro h, one has first 
to solve Eq. (J33|) to get Xond, then Eqs. (JMJ) and (J36J) to 
obtain x*" and x* nd horn which A7V^ UC and AN* nd are 
deduced by using Eq. (|32]l. From the value of x*, one 
can also directly evaluate the two slow-roll parameters 
ei* and 62*. Let us notice that some of the expressions 
above can be ill defined if p = 2. In this case, Eqs. 
and (|3l>|) should be rederived from the start and one can 
show that it then always leads to well-defined expressions. 
The rest is the same as for the large field models (see 
Sec. Ill Cj ). proh being in one-to-one correspondence with 
the value of AA^, through Eqs. (fTU)) and fllUl) once a value 
of w Te h has been chosen. It is worth emphasizing again 
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FIG. 2: Reheating consistent slow-roll predictions for the large field models in the plane (£1,62)- The two blue dot-dashed 
contours are the one- and two-sigma WMAP3 confidence intervals (marginalized over second order slow-roll) while the pink 
solid contours are the one- and two-sigma WMAP7 ones. As in Fig. [TJ the annotations trace the energy scale at which the 
large field reheating ends and correspond to log(<jy 4 T re h/GeV). The solid line represents the model p > 1. This confirms that 
there now exists a lower bound on the value of <?* T re h ( see Sec. IIII|) . 



that, in order to derive these results, no assumption has 
been made about the reheating epoch which is entirely 
characterized by p rc h and w) re h- 

In Fig. [31 we have represented the slow-roll predictions 
for an assumed generic value of p = 3 and various val- 
ues of w Te h ranging from —0.2 to 0.8. The annotations 
are the values of log(gy 4 T rc h/GeV) while the color scale 
traces the values of fi/M Pl . For small field models, the 
reheating energy scale is all the more so constrained that 
w} re h and /j. are small. In fact, these plots show that w Te h 
and n are degenerated: it is possible to render compatible 
a low value of w rc h provided (i is super-Planckian. Con- 
versely, small values of /j, can be made compatible with 
the data for a high energy scale reheating if w rc h < 1/3, 
or a low energy scale reheating for w re h > 1/3. In the 
next section, we perform a full analysis of the small field 
models, reheating included, in view of the WMAP7 data 
to quantify the above claims in terms of posterior prob- 
ability distributions. 



III. INFERRING REHEATING FROM CMB 
DATA 

In view of the previous results, the correct way to dis- 
cuss how well the CMB data constrain a set of known in- 
flationary models is to perform a Bayesian analysis of the 
data given the model parameters, including the reheat- 
ing. Notice that this is different than constraining the 
slow-roll parameters, or the spectral index and tensor- 
to-scalar ratio, which only encode the shape of the pri- 
mordial power spectra and know nothing about reheating 
(whereas a model of inflation does). 



A. Exact numerical integration 

The numerical exact integration method has been in- 
troduced in Refs. [13, S3, [H, [53[ and consists of the 
computation of the primordial power spectra assuming 
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FIG. 3: Reheating consistent slow-roll predictions for the small field models with an assumed generic value p = 3 and for 
w la h — —0.2 (top), Wreh = (middle) and w re h = 0.8 (bottom). The right panels display the corresponding predictions in 
the plane (n s ,r). The reheating has a strong influence for low values of both w re h and n/M p[ . As in Figs. [T] and the 
annotations give the values of log(<7y 4 T re h/GeV). It is also interesting to notice that the sequence of successive T re h is switched 
for w r eh = 0.8 (in fact for w la h > 1/3), i.e. large reheating temperatures correspond to smaller spectral indices for w Ic h — 0.8 
(i.e. w rc h > 1/3) while, for w rc h = or w rc h = —0.2 (i.e. w rc h < 1/3), they correspond to larger n s . 



only General Relativity and linear perturbation theory. 
Therefore, the only model parameters are the ones ap- 
pearing in the inflaton potential together with the re- 
heating parameter i? ra d, for the very reasons explained 
in Sec. HU The numerical integration of the inflation- 
ary perturbations sets up the initial conditions for the 



subsequent cosmological perturbations from which the 
CMB anisotropics arc deduced. For this purpose, we have 
used a modified version of the CAMB code [13] coupled 
to a Monte-Carlo-Markov-Chain (MCMC) exploration 
of the parameter space implemented in the C0SMDMC 
code Hi and given the WMAP7 data 0-1111 . Con- 
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cerning the standard cosmological model, we have as- 
sumed a flat ACDM model having five parameters: the 
density parameter of baryons, f2b, of cold dark matter 
fidm, the Hubble parameter today Hq, the optical depth 
t encoding the redshift at which the Universe reion- 
ized, and the nuisance parameter A sz encoding the rel- 
ative amplitude of the diffuse Sunyaev-Zel'dovich (SZ ) 
effect compared to the anal ytic al model of Ref. [56j . 
In fact, as discussed in Ref. [55[, it is more convenient 
to sample the cosmological parameter space along the 
rescaled quantity (Q,^h 2 , Q^ m h 2 , r, 9, A sz ) where H = 
100/ikm/s/Mpc and 6 measures the ratio of the sound 
horizon at last scatte ring to the angular diameter dis- 
tance. Following Ref. [4a], we have included the lensing 
corrections on the temperature and polarization power 
spectra, and, to limit parameter degeneracies, completed 
the WMAP7 data with the latest Hubble Space Tele- 
scope (HST) bound on Hq [53]. Concerning the primor- 
dial parameters, they are now provided by our inflation- 
ary model parameters, up to some observationally con- 
venient rescaling. For instance, we will prefer to sample 
on P*, the amplitude of the scalar perturbation at the 
pivot scale, rather than on the potential normalization 
M, both being in one-to-one correspondence. Similarly, 
it is more convenient to sample the reheating era over the 
parameter P, 



R = R 



rad 



1/4 
lend 



(38) 



rather than P ra d- As can be seen by comparing Eq. ([15]) 
and the following exact expression: 

\*R = N T -N*+N Q + \\*U^±^-\, (39) 

contrary to R, the values of P ra d explicitly depend on 
P*. This would induce unwanted correlations between 
-R ra d and P* which are therefore avoided by sampling the 
reheating over R. Notice that since R and P ra d differ 
by a factor p cn d, they are also in one-to-one correspon- 
dence once the model of inflation is specified [27| . In 
order to perform the MCMC analysis, we still have to 
specify the prior probability distributions. Concerning 
the cosmological parameters, we have chosen wide flat 
priors around the preferred posterior values obtained by 
the WMAP team [45|. The reheating energy scale being 
unknown, we assume a flat prior on InP whose exten- 
sion is given by the consistency conditions mentioned in 
Sec. Ill Bl Reheating should occur before nucleosynthesis 
and after the end of inflation while the positivity energy 
conditions imply — 1/3 < u> rc h < 1- As a result, we take 
a flat prior for In R in the range [27 



In 



'1/4* 
pnuc 

HI 



< In R < — In 
3 




3 M„ 



The lower bound is approximately ~ —47 for p nuc = 
lOMeV whereas the upper bound depends on p cn d, and 



thus on the other inflationary model parameters. Fi- 
nally, we have chosen a flat prior on the logarithm of P* 
around the value giving the right amplitude of the CMB 
anisotropies: 2.7 < In (l0 lo P*) < 4.0. Let us notice that 
since P* is well constrained, this translates into an upper 
bound on ln(p^/M p J < —5.5 (analogous to the upper 
bound on in the slow-roll approximation) that will 
therefore be inherited by In R so that the maximal value 
of the upper bound is ~ 8.3. The other prior choices 
on the primordial parameters are those concerning the 
inflaton potential and will be specified later. 

In the following, we perform the WMAP7 data anal- 
ysis along those lines for both the large and small field 
models. As the first step, we sample over the rescaled re- 
heating parameter In R without any assumptions on w re h- 
We show that it is actually constrained for all models. 
As can be checked in Eq. (ITUl) . it means that the CMB 
data restrict the a priori possible values of AN and w re h- 
Conversely, this result shows that not including the re- 
heating parameter when constraining inflationary mod- 
els is no longer a reasonable option. For the second step, 
we break the degeneracy between w re h and AA~ and as- 
sume that w re h takes its natural values for the large field 
models given in Eq. (f2"3|) . or choose a specific value in 
the small field models. These reasonable extra assump- 
tions translate the bounds on InP into a lower limit on 
/9 ron and/or gl^T^. Unless specified, we have stopped 
the MCMC exploration according to the R-statistics [58[ 
implemented in COSMOMC such that the difference in vari- 
ances between the different Markov chains does not ex- 
ceed a few percent. Typically, this corresponds to a set of 
300 000 to 500 000 samples depending on the underlying 
model of inflation. 



B. Large field models 

The potential for the large field models is given in 
Eq. (|2"Tj) . Together with P* and In R, there is only one ad- 
ditional primordial parameter p for which we have chosen 
a flat prior in the range p € [0.2, 5]. The upper bound is 
motivated by the previous constraints on large fields 27] 
whereas the lower one is a theoretical prejudice associated 
with the non-naturalness of extremely small values of p in 
any field theory. The marginalized posterior distributions 
for the sampled and derived cosmological parameters are 
represented in Fig. [4] for the two prior assumptions de- 
tailed in the following. The solid lines are without any 
assumption on the reheating whereas the dashed ones are 
under the natural equation of state iD ro h = (p— 2) / (p+2). 
These probabilities are compatible with the one already 
derived in the literature [45|, [59| up to slight shifts coming 
from changing the reheating assumptions. This is the re- 
sult of some tension between the large field models which 
generically predict a large tensor-to-scalar ratio r and its 
non-observation. Being more restrictive on the reheating 
gives less flexibility to the model such that Ho and tth are 
slightly shifted to compensate for the too high r values. 
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FIG. 4: Marginalized posterior probability distributions for 
the base and derived cosmological parameters in large field 
inflation. The black solid lines are without any assumptions 
on the large field reheating whereas the red dashed ones are 
under the prior w Ic h = (p — 2)/(p + 2). They have been 
represented only when they differ with respect to the former. 
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FIG. 5: Marginalized posterior probability distributions (solid 
lines) and mean likelihoods (dotted) for the large field model 
primordial parameters. This is without assumption on the 
reheating era. Notice the lower bound on the reheating pa- 
rameter In R which correlates with the potential power p (see 
also Fig. [6}. The energy scale at the end of large field inflation 
is also constrained. 




0.5 


1 


1.5 


2 2.5 
P 


3 


3.5 


4 




-27 


-26 


-25 


-24 , -23 
ln(p /M ,) 
vl end pi' 




-22 


-21 



FIG. 6: Two-dimensional marginalized posterior probability 
distribution (point density) in the plane (p, In R) and its one- 
and two-sigma confidence intervals. Correlations with the en- 
ergy scale of large field inflation are traced by the color scale. 



In Fig. we have plotted the marginalized proba- 
bility distribution for the large field primordial param- 
eters without assumption on the reheating. It is partic- 
ularly interesting to compare these plots to Fig. 18 of 
Ref. j27[ since this allows us to see the improvements on 
the parameter constraints coming from the passage from 
WMAP3 data to WMAP7. In addition to the expected 
constraints on P*, we find the 95% confidence limit 



p<2.2, 



(41) 



suggesting that 4> 2 inflation may now be considered un- 
der pressure. Let us emphasize that this result is robust 
against any possible reheating evolution since marginal- 
ized over In R. Concerning this last parameter, we find a 
95% lower bound: 



lni? > -28.9. 



(42) 



In fact, as can be checked in Fig. [HI these two parameters 
are correlated together and also with p C nd- These corre- 
lations can be understood as follows. From Eq. (|39)) . the 
quantities lni?, p and In (p cn a/M^^ are related by the 
formula 



In 



( !>■ 



'end 



= In (l28vr 2 F») - 2N + 2 In R - 



pi + n s 
~2 1-n = 



In 



8p(l-n s ) 



P- 



(43) 



where P, and n s are well constrained quantities. As a 
result, at fixed lni?, the larger the p values, the lower the 
energy scale at the end of inflation has to be, which is 
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FIG. 7: Marginalized posterior probability distributions (solid 
lines) and mean likelihood (dotted lines) for the large field pa- 
rameters when it) ro h = (p — 2)/(p + 2). This extra-assumption 
on reheating yields to tighter constraints than in Fig. [5] In 
particular, we find p rc h > 17.3 TeV at 95% confidence level, 
as well as p < 2.1. 



exactly what is observed in Fig. [6j Of course, p cannot 
be too large since, in this case, the tensor-to-scalar ratio 
r increases and rapidly becomes incompatible with the 
CMB data. In Fig. [51 we also observe that the smaller 
p, the larger the allowed range of variation of lni?. The 
upper limit on In i? does not depend on p and just comes 
from the upper limit on the energy scale of inflation 
(lni? < 8.3). On the other hand, the lower limit strongly 
depends on p and represents a non trivial result. This 
expresses the fact that, for a given p, there are values 
of In i? for which there is no way to obtain, at the same 
time, a consistent reheating epoch and CMB predictions 
compatible with the data. From this effect, we also get 
the energy scale of large field inflation, and at two-sigma 
level 

4.4 x 10 15 GeV < pj^ < 1.2 x 10 16 GeV . (44) 

The upper limit just comes from the constrain on the 
energy scale of inflation while the lower limit originates 
from the fact that p cannot be too large (recall low values 
of Pend means large value of p). 

Assuming now that the reheating proceeds according 
to Eq. (|2"3")l , one obtains the marginalized posteriors plot- 
ted in Fig. [7] and Fig. [5] To be consistent, we have mod- 
ified our prior on p by assuming a flat distribution in 
p g]1,5]; the case p = 1 is a limiting case that may 
be problematic. Indeed, values of p < 1 would induce 
w ro h < — 1/3 and inflation would not stop. Moreover, 
instead of using lni?, we have used Eqs. ([TBI and (l38l) to 
sample the parameter space over ln(pteh/M^ ) and from 
a flat prior in [ln(p nuc /M^),ln(p end /-W^)]. The upper 
bound on the p posterior is slightly tighter due to the re- 
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FIG. 8: One- and two-sigma marginalized limits in the plane 
[p, ln(/9 re h/.Mp ,)] for the large field models with w Te h = (p ~ 
2)/(p + 2). The point density traces the associated two- 
dimensional posterior while the color map shows correlations 
with the energy scale at which inflation ends. 



striction made over the reheating: we find the two-sigma 
limit p < 2.1. For the same reasons, the energy scale 
of large field inflation is a bit more constrained, and the 
two-sigma range becomes 

5.2 x 10 15 GeV < pj^ < 9.1 x 10 15 GeV. (45) 

Certainly, the more interesting result is the lower bound 
on the reheating energy scale. At 95% of the confidence 
limit 

pt£> 17.3 TeV. (46) 

The correlations between these three parameters are rep- 
resented in Fig. [8] and have the same origin as the ones 
displayed in Fig. [SI up to the change of variable i? to p ro h. 
In particular, we see that, at a fixed value of p, the con- 
straints on /9 rc h are tighter for p < 1.5 than for p ~ 1.5. 
This comes from the fact that the reheating is well con- 
strained for a negative mean equation of state, which 
precisely corresponds to p < 2 [see Eq. (1231) ]. The change 
of behavior around p — 1.5 comes from this effect com- 
bined with a two high tensor-to-scalar ratio when p > 2. 
Let us also emphasize that Eq. (1461) is marginalized over 
all large field models. Coming with a theoretical prefer- 
ence for a given value of p can lead to stronger bounds, 
as for instance if p > 1 or p = 2 (see Fig. [8]). Finally, 
one can check that the bounds found in this section are 
compatible with the expectations we have derived from 
the slow-roll predictions of Sec. Ill CI 
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FIG. 9: Marginalized posterior probability distributions (solid 
lines) and mean likelihoods (dotted) for the small field model 
primordial parameters. This is without assumption on the 
reheating era. Notice again the lower bound on the reheat- 
ing parameter lni? and the slightly favored super-Planckian 
values of fi. Correlations are displayed in Fig. I10I 



In the next section, we perform a similar analysis for 
the small field models. 



C. Small field models 

The small field model potential of Eq. ([31]) involves 
an extra parameter compared to large fields which is the 
VEV (i. The scale of this parameter being unknown, 
we have chosen a flat prior on log(///Af pi ) in the range 
[—1,2]. With the lower and upper limits being chosen 
only for numerical convenience, one should keep in mind 
that the physical values of /i may be larger or smaller. 
The important point is however that such a prior includes 
both sub-Planckian and super-Planckian values without 
prejudice. Concerning the potential index p, we have 
chosen a flat prior in the range p C [2.4,10]. The up- 
per bound is arbitrary whereas the lower bound excludes 
p = 2 since this model is a special case [27]]. In fact, there 
exists another reason that justifies the above choices. In 
the limit n/M Fl ^> 1, one can show, using a perturbativc 
expansion in M pl //j,, that the two first horizon flow func- 
tions ei and ei become independent from /i andp, namely 
€i = (4AiV» + 1) _1 and e 2 = 4ci. Therefore, it would be 
useless to take a larger upper bound on the fi prior since 
the corresponding physical predictions are no longer af- 
fected by this choice. All the other priors, both on the 
cosmological and primordial parameters, have been cho- 
sen as for the large field exploration. For the sake of 
clarity, we have not represented the marginalized poste- 
riors for the cosmological parameters. Contrary to the 
large field models, these posteriors are almost the same 



whatever we assume for the reheating. The reason is that 
small field models do not have a tendency to produce a 
high tensor-to-scalar ratio. There is therefore no need 
for the cosmological parameters to compensate for such 
an effect and they decouple from the details of the infla- 
tionary and reheating phases. Finally, the cosmological 
parameter posteriors in small field inflation end up being 
very similar to the dashed curves (or solid curves when 
they are absent) plotted in Fig. O 

In Fig. IH1 and Fig. [TU1 we have plotted the marginal- 
ized posterior probability distributions (one- and two- 
dimensional, respectively) for the primordial small field 
parameters without assumption on the reheating. Again, 
we find the WMAP data to give a lower limit on the re- 
heating parameter, at two-sigma level 

lni? > -23.1. (47) 

Concerning the parameters \i and p, they are not con- 
strained. However, the posteriors of Fig. [§] clearly show 
a tendency to favor super-Planckian values of /x together 
with large values of p. As can be seen in Fig. [TUl since /j, 
can take arbitrarily low values, the energy scale of small 
field inflation is not constrained from below. We find 
only the consistency condition that < 9 x 10 15 GeV 
from the P* limits. The correlations between fj, and p can 
be understood from Sec. Ill Dl and come from the require- 
ments of having the right spectral index. The bound of 
Eq. (|4T|) has the same origin but through the selection 
of the favored AiV* values. More details on these effects 
can be found in Refs. [13, HzJ. 

As we did for large fields, we now assume an equation 
of state parameter for the small field reheating. Con- 
trary to large field models, we no longer have a relation- 
ship between u> re h and the inflationary potential and one 
may only assume fiducial values ranging from —1/3 to 
1. Again, the MCMC has been now sampled directly 
over ln(/O re h/Afp, ) rather than over In R by making use of 
Eqs. (fT3"j) and l[38"]). The resulting one-dimensional poste- 
riors for the primordial parameters have been plotted in 
Fig. QT] for the four values w T eh = —0.3, —0.2, —0.1 and 
w lc h = 0, and only when the posteriors are affected by 
this choice. Comparing Figs. |H] and Fig. [TT] shows that 
the posteriors of /it, p (and P*) are mostly independent of 
the details of the reheating. The only quantity changing 
accordingly is p rc h ■ This is not surprising since at a given 
R, changing the values of w Ic h modifies the number of c- 
folds the Universe reheated. As a result, the constraint 
on R for small fields translates into a lower bound on /o re h 
but only when w} re h is small. At 95% of confidence, we 
find 

w ieh = -0.3 => p l J2 > 8.9 x 10 5 GeV, 

(48) 

w rch = -0.2 => pV* > 3.9 x 10 2 GeV, 

while higher values of w Te h do not constrain p re h more 
than the prior p rc h > p n uc- Physically, these results 
can be fully understood from Fig. [3J In particular, the 
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FIG. 10: One- and two-sigma contours of the two-dimensional marginalized probability distributions for small field inflation. 
The point density traces the associated two-dimensional posterior while the color scale shows correlations with third parameters. 
Notice the correlation between fi and p, as well as with In R. 




FIG. 11: Marginalized posterior probability distributions 
(solid lines) and mean likelihoods (dotted lines) for the small 
field model. These posteriors assume a constant equation of 
state with (from bottom to top) w Ie h = —0.3 (black lines), 
—0.2 (red lines), —0.1 (blue lines) and w Ic h = (green lines). 
Only p re h is affected by such prior choices (bottom right 
panel). The energy scale of reheating is all the more con- 
strained from below than w Ic h is small. 
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FIG. 12: One and two-sigma marginalized limits in the plane 
[log(/i/M P! ),ln(p reh /Mp i )] for the large field models with 
w rc h = —0.3. The point density traces the associated two- 
dimensional posterior while the color map shows correlations 
with the energy scale at which inflation ends. 



fact that constraints on p vc h can be derived for nega- 
tive values of u> r ch is apparent from these plots. As ex- 
pected, correlations between /9 rc h and the other param- 
eters are similar to the ones associated with \nR. In 
Fig. I12[ we have represented the two-dimensional pos- 
terior and its one- and two-sigma contours in the plane 



[log(/Lt/M P! ),ln( / o re h/M p l )]. The color scale traces corre- 
lations with the energy scale of small field inflation. We 
recover the slightly disfavored values of sub-Planckian 
vacuum expectation values. They clearly remain accept- 
able but only if the reheating ends at a high energy. 
Again, the previous results are compatible and easily un- 
derstandable with the slow-roll predictions of Sec. Ill Dt 
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To conclude this section, we have found that the cur- 
rent WMAP data give non-trivial information on the re- 
heating era in the two most considered classes of proto- 
typical inflationary models (large and small fields). Con- 
versely, when the goal is to use CMB data to constrain 
these models, including the reheating into the marginal- 
ization is definitely no longer an option. 



D. Importance sampling from slow-roll bounds 



The formulas derived in Sec. Ill Bl link the reheating 
and inflationary model parameters to the spectral index 
and tensor-to-scalar ratio. As a result, they could be 
used to extract constraints on the reheating energy scale 
from some already derived constraints on the slow-roll 
parameters by using importance sampling [551 ] . However, 
this method is highly inefficient since most of the favored 
(n s , r) values do not necessarily correspond to a consis- 
tent reheating model in a given inflationary framework. 
However, it clearly illustrates how knowledge on the pri- 
mordial power spectra shape, complete with an inflaton 
potential, translates into some information on the en- 
ergy scale at which the reheating ends. For this reason, 
we briefly discuss this method in the following although 
we prefer an exact numerical integration as performed 
above. 

Assuming an inflationary model, with a known poten- 
tial V(<j>), one can generate any spectral index and tensor- 
to-scalar ratio such that, at leading order in slow-roll, 



r = 16ei», n s — 1 = — 2ei* — 62*, (49) 



with 



1 (VI 



£2* — — 2 



Yl_(K 



(50) 



A prime here is understood as a derivative with respect to 
the field <fi. The star still refers to the time at which the 
scale under consideration crossed out the Hubble radius 
during inflation. Inverting Eq. (|4*91 , or Eq. (I5T)]) , gives the 
value of 4>* (and eventually other parameters) leading to 
the required couple (n s ,r). For instance, in the large 
field models, one would find 



'Sen 



1 pi V 2* 



p = 4 — 

£2* 



(51) 



From </>* one can derive AiV, from the slow-roll trajectory 
and the value at which inflation stops, (f> cn d- For instance, 
for large field inflation, one would obtain 



AA* 



1 - £l* 
£2* 



(52) 



Eqs. (50) and (fT3j): 

3 + 3iD roll 



In 



' 1/4- 

Proh 



1 - 3w re h 

1 + 3w roh 




(53) 



e lend V* J 



provided w TC h 7^ 1/3. As expected, if the reheating is 
radiation dominated, it cannot be distinguished from the 
usual radiation era and p re h cannot be inferred. Clearly, 
for a given set (ei*, C2*, P*), in an assumed model of 
inflation, the right-hand side of Eq. (|53|) is uniquely 
determined and hence is /? r ch- As already mentioned, 
such a method is not well suited: picking up a ran- 
dom (ei*,£2*,P*) compatible with the power spectra 
shapes usually predicts a value of p rc \ l which is either 
incompatible with BBN, or with the underlying model, 
i.e. p ro h > Pend- The reason is that inflationary physics 
is much more than Taylor expanding a potential and fit- 
ting the power spectra shape. In order to solve this issue, 
the way out is to perform an exact numerical integration 
of the inflationary perturbations, including the reheating 
parameter, as we previously did. This method has also 
the advantage to free ourselves from any assumption on 
the equation of state parameter w re h- 



IV. CONCLUSION 

We now conclude our investigation by revisiting our 
main results. The most important conclusion is that, 
both for large and small field scenarios, the reheating 
parameter lni? is now constrained by CMB data. The 
physical origin of this result is clear. For fixed physical 
length scales today, a change in In R modifies the location 
of the CMB observable window along the inflationary po- 
tential, which is possible only for a limited range of In R 
given the data accuracy. This conclusion is general and 
does not depend on the details of the reheating epoch. 
However, if one assumes a model for reheating, typically 
if one chooses a specific value of the mean equation of 
state, then it becomes possible to express the constraints 
mentioned above as limits on the energy density at the 
end of reheating or, equivalently, as constraints on the re- 
heating temperatures. This leads to Eqs. P5t . and (|4"5j) . 
These results are of particular interest for the supersym- 
metric extension of either large or small field inflationary 
models. Indeed, our result limits the reheating tempera- 
ture from below whereas gravitinos production gives an 
upper bound [32h42| : T rc h < 10 4 TeV, where T re h is given 
in Eq. (f29|) . Therefore, assuming ~ 200 in the large 
field case, we now have an allowed range of variation for 
the reheating temperature given by 



The energy scale at which reheating ends stems from 



6 TeV < T rch < 10 4 TeV. 



(54) 
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By including the reheating parameter in our analy- 
sis, we can marginalize over all reheating history to infer 
the inflationary parameter values in a robust way. For 
large held scenarios, we hnd that the power index is up- 
per limited by p < 2.2, at 95% confidence limit. This 
means that the prototypical model of inflation, namely, 
massive chaotic inflation, is now under pressure. Simi- 
larly, the small field models with sub-Planckian vacuum 
expectation values /j, are slightly disfavored. In fact, [i is 
correlated with the potential power p, as represented in 
Fig. [TUJ Without marginalization over p, large values of 
p > 6 are actually necessary to allow the sub-Planckian 
values of n to be inside the 95% contour. If one has a the- 
oretical prejudice for /i < M rl , and for reasonable values 
of p < 6, then small field models can also be considered 
under pressure. 

Since the constraints on the reheating parameters are 
directly related to the ability of the data to determine 
the observable parameters, one could consider more data 
sets than the WMAP7 data. In fact, since solely the ac- 
curacy on the primordial parameters matters, data sets 
improving the constraints on the standard cosmological 
parameters do not change the reheating bounds. On the 
other hand, small scale CMB experiments may be deci- 
sive but, as mentioned in Ref. [45|], they do not give a 
significant improvement on the determination of n s and 
r due to their low accuracy at large multipoles. We have 
indeed tested that our limits do not change by including 
the Baryonic Acoustic Oscillation [6(3] or the Arcminute 
Cosmology Bolometer Array Receiver data [6l| in our 
analysis. On the other hand, since the future Planck data 
are expected to improve the bounds onn s , r and even on 
new primordial observables, we should get unprecedented 
information on the inflationary reheating era. 
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Appendix A: Reheating from inflaton decay 

For large field models, we have used the fact that 
Wroh = (p — 2)/(p + 2), a well-known formula established 
for the first time in Ref. [l7j- However, this result as- 
sumes that the inflaton field is not coupled to other fields, 
a hypothesis that, if acceptable at the early stages of the 
reheating phase, cannot be maintained if one is inter- 
ested in the transition to the radiation dominated era. 
The purpose of this Appendix is to take into account the 
coupling of the inflaton field with the rest of the world 
and to study its impact on the value of u) rc h. As we 
show in the following, when the coupling is considered, 
the value of w> re h does not significantly deviate from the 
equation given above. 




E-fold number N 




E-fold number N 

FIG. 13: Top panel: Evolution of the inflaton (dotted green 
line) and radiation (solid red line) energy densities during the 
reheating epoch as a function of the number of e-folds since 
the beginning of large field inflation with p = 2 (N T ~ 61). 
The inflaton decay rate has been chosen to be F ~ 1.375 x 
10 9 GeV corresponding to a reheating temperature of T rc h — 



3.2 x W 13 g 



-1/4, 



GeV. The total number of e-folds during the 
reheating epoch is iVreh = 5. Bottom panel: the instantaneous 
equation of state (black solid line), w rcn = (-Pinf+p7/3)/(pmf+ 
p 7 ), during the reheating era and going to 1/3. The dashed 
red line represents w Ic h- Its value when the instantaneous 
equation of state has reached 1/3 is ~ 0.08 so that, even when 
the number of e-folds during reheating is small, deviations 
from u) re h = never exceed 8%. 



A simple and standard way to model the physical sit- 
uation where the inflaton field decays into radiation is to 
write the Klein- Gordon equation as 17] 



dV 

(3H + r)<t>+ — = 0, 



(Al) 



with r being the inflaton decay rate. The evolution equa- 
tion for the radiation energy density is modified accord- 
ingly in order to ensure that the total energy density is 
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tog(r) 



-0.2 -0.1 




FIG. 14: Marginalized posterior probability distributions 
(solid lines) and mean likelihoods (dotted lines) for the three 
first slow-roll parameters and for the energy scale of inflation 
(k=1/M p1 ). 



conserved. One obtains 



2 P r 

7T2H 



Pint, 



(A2) 



where p- ln f = <f> 2 /2 + V(<f>). Of course, these two formu- 
las must be supplemented with the Friedmann-Lemaitre 
equation, H 2 = (/Qj n f + Pj) /(3 M^1 . We have numeri- 
cally integrated Eqs. (|A1[) and (|A2|) for large field infla- 
tion with p = 2 and the corresponding evolutions of p m f , 
Prad, w re h and w} re h are displayed in Fig. [T3l 

There are some subtle issues if one wants to numeri- 
cally evaluate the value of w re h- For instance, the time 
at which one considers the Universe to be reheated is not 
very well defined since we have a smooth transition, and 
this affects the precise numerical determination of w le h- 
Indeed, one could consider that reheating is completed 
when t ~ r _1 or when w rc h — 1/3 ± <5 for a given 6. 
These choices (for instance the precise value of S) lead to 
different values of the mean equation of state. Moreover, 
the numerical calculation itself can be difficult since one 
has to integrate a rapidly oscillating function. All in all, 
we find that, for A^ rc h ~ 5 (as in Fig. Q2]), w Te h < 0.08, 
that is to say a value close to 0. Let us also notice that 
when the value of -/V ren increases (i.e. when the reheat- 
ing temperature decreases), one expects this value to be 
even less than the limit quoted before (but this regime is 
numerically difficult to follow since the equation of state 
rapidly oscillates during a long time). Therefore, we con- 
clude that Eq. ([23")) can reasonably be trusted and this 
justifies the approach used in this paper. 



FIG. 15: Marginalized posterior probability distributions 
(solid lines) for the spectral index and the tensor-to-scalar 
ratio in the slow-roll approximation. The primordial power 
spectra are expanded at second order in slow-roll. 



Appendix B: Slow-Roll posterior distributions 

In Sec. |nl we used the one- and two-sigma contours in 
the planes (n s , r) and (ei, e 2 ) obtained from the WMAP7 
data. In this section, for the sake of completeness, we also 
provide the resulting one-dimensional marginalized pos- 
terior probability distributions on those four parameters 
together with the third slow-roll parameters €3 and the 
energy scale of inflation. Let us notice that our analysis, 
even for n s and r, assumes the second order slow-roll ex- 
panded primordial power spectra for both the scalar and 
tensor perturbations. 

In Fig. 1141 we present the marginalized posterior prob- 
ability distributions for the three first slow-roll param- 
eters and for the Hubble parameter during inflation 
(H[ n f — For the first slow- roll parameter, we have 

assumed a Jeffreys' prior in the range [10 _5 ,0.2], i.e. a 
flat prior on log(ei) in the range [—5,-0.7], as appro- 
priate for a parameter the order of magnitude of which 
is unknown. For the two next slow-roll parameters, we 
have chosen flat priors in [—0.2,0.2]. The upper limit on 
ei directly comes from the level of primordial gravita- 
tional waves. The parameter e 2 is well constrained while 
63 remains unbounded. There exists an upper limit on 
the energy scale of inflation which, as for the first slow- 
roll parameter, directly comes from the non-observation 
of the primordial gravitational waves. In Fig. 1151 the 
spectral index and the tensor-to-scalar ratio have been 
plotted and obtained by importance sampling from the 
slow-roll constraints. Indeed, at second order in slow-roll, 
one has 

n s = 1 - 2ci - e 2 - 2ej - (2C + 3) ei e 2 - Ce 2 e 3 

_2 



r = 16ei 



1 + Ce 2 + ( C - — + 5 ) £l e 2 



(Bl) 



where C is a numerical constant, C 



-0.7296. As can 



be seen on these plots, the spectral index value n s = 1 
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is disfavored but not excluded when marginalizing over limit on the tensor-to-scalar ratio r. 
second order slow- roll. As for ei, there is only an upper 
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